Boson Mott insulators at finite temperatures 
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We discuss the finite temperature properties of ultracold bosons in optical lattices in the presence 
of an additional, smoothly varying potential, as in current experiments. Three regimes emerge in 
the phase diagram: a low-temperature Mott regime similar to the zero-temperature quantum phase, 
an intermediate regime where MI features persist, but where superfluidity is absent, and a thermal 
regime where features of the Mott insulator state have disappeared. We obtain the thermodynamic 
functions of the Mott phase in the latter cases. The results are used to estimate the temperatures 
achieved by adiabatic loading in current experiments. We point out the crucial role of the trapping 
potential in determining the final temperature, and suggest a scheme for further cooling by adiabatic 
decompression. 

PACS numbers: 03.75.Lm,03.75.Hh,03.75.Gg 



Ultracold gases in optical lattices are currently a major 
topic in the field of ultracold atoms, both experimentally 
and theoretically (see the reviews P, [3j S ^^^^ 
experiments are mostly carried out with Bose gases, and 
interpreted in terms of the zero-teniperature phase di- 
agram of the Bose- Hubbard model [J, H, Q- However, 
actual experiments inevitably take place at finite tem- 
peratures, which leads, e.g., to residual number fluctua- 
tions. In view of its importance for applications such as 
the controlled generation of entangled states or the 
study of quantum magnetism |3|, this issue has rece ntly 
received increasing interest [1, 0, 0, S, S [13, [O, [13, [l3| . 
In this paper, we first discuss how the phase diagram 
of a ultracold Bose gas in an optical lattice is modified 
at finite temperatures. We identify a melting tempera- 
ture T* « 0.2 U/k-Q, above which the system is entirely 
thermal, and a much lower critical temperature Tc ~ zJ, 
above Mott regions survive, but superfluidity is absent. 
Here z is the coordination number, U the on-site inter- 
action energy, and J the tunneling energy. We derive 
analytically the thermodynamics of the Mott phase for 
T < T* , including particle-hole corrections. Finally, we 
model the procedure adopted in current experiments to 
produce quantum gases in optical lattices, where a Bose- 
Einstcin condensate is slowly loaded into the optical po- 
tential in the presence of an "external" , smoothly varying 
trap. Assuming the loading is done adiabatically, we es- 
timate the final temperatures that can be reached with 
this procedure. We point out the crucial role of the exter- 
nal potential in determining the final temperature, and 
suggest a scheme for further cooling by adiabatic decom- 
pression. 

For future reference, we first recall the phase diagram 
for zero temperature. For J = 0, the ground state is a 
product of identical Fock states, with uq atoms at each 
site depending on the chemical chemical potential fi ac- 
cording to no = Integer(///C/) + 1. When a weak tun- 
neling between nearest neighbors is allowed, the ground 



3 
2.5 
2 

i,.5 

1 

0.5 






Superfluid 


■ ' Normar^*''''^ 




. Thermal 


QuantuTTT 


Insulator 


Insulator ^ 







0.05 



0.1 
zJ/U 



0.15 



0.2 



FIG. 1: Phase diagram for an ultracold Bose gas in an opti- 
cal lattice at finite temperatures (ksT — U/20). The finite- 
temperature Mott lobes boundaries are shown as solid lines, 
and the zero-temperature boundaries as dotted lines for com- 
parison. With increasing temperature, the superfiuid regions 
in between the Mott lobes shrink in width and shift to larger 
tunneling energies. For temperatures zj < T, the SF-MI 
transition is replaced by a smooth crossover between two Mott 
domains across a normal phase. This generic phase diagram 
persists up to the melting temperature T* ~ 0.2 U where the 
Mott domains vanish. 



state remains a Mott insulator with suppressed num- 
ber fluctuations but only in well-defined domains ( "Mott 
lobes") in the J — fi plane (dashed lines in Fig.[T|). Out- 
side these domains, the MI is unstable towards a delo- 
calized, superfiuid state (SF), that connects to a conden- 
sate in the Bloch state q = when J ^ U . In current 
experiments, an external harmonic potential of the form 
Vt{t) = ^mcu^r^ is usually superimposed on the lattice 
potential. As a result, the density profile is in general 
not uniform everywhere, but has the shape of a "wed- 
ding cake" formed by successive layers ("Mott shells") 
with integer densities [13]. If this potential varies slowly 
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over a few lattice sites, the local density approximation 
(LDA) can be employed. It obtains the (coarse-grained) 
number nij[/j,, T] and entropy Sh[/x, T] densities calculated 
as in the homogeneous case, but for a local chemical po- 
tential fiQ — VT(r) = ^ioc(i")- In the following, we first 
compute the phase diagram for uniform systems, and use 
the LDA to describe trapped gases. 

Within a Mott domain with hq bosons per site, the 
lowest-lying excited states are "particle" and "hole" 
states, where a supplementary particle is added (respec- 
tively removed) to the background value uq with a free 
energy cost Uhq — fi (rcsp. ^ ji — U {uq — 1)). Other occu- 
pation numbers correspond to free energies at least of or- 
der U, and are therefore suppressed by a thermal factor ~ 
g-/3(7 ^ ^ This allows to use what we call the "particle- 
hole approximation" (PHA) [1, i, [H, [H, [I^, El , which 
consists in truncating the on-site Hilbert space to the 
states no, uq — 1 and tiq + 1. As we will see, this ap- 
proximation, valid for temperatures fc^T U , allows 
to describe quite extensively the physics of the insulating 
phase at finite temperatures. 

As a starting point, we first neglect tunneling be- 
tween the wells completely and consider the thermo- 
dynamics of a uniform system of isolated wells, where 
the atoms interact locally through the on-site interac- 
tion. Since the wells are independent, the global parti- 
tion function factorizes into a product of identical on-site 
partition functions zq = X^n^'^P ~ Z^")]? where 

E{n) = Un{n - l)/2 and /? = l/fceT. Keeping only the 



particle/hole excitations with energies 

^o'^^'' ~ U{ni) ^1)7 the mean density ai 
density fluctuations take the simple forms 



(qp) 



U no and 



no « 710 + - S(-)) /zo, (1) 

Var(n)o « (i?(+) + i?(-)) /z^, (2) 

where = e''['^-'^"«l and B^"' = g/3[c/(no-i)-M] are 

Boltzmann factors corresponding respectively to adding 
or removing a particle from the "background" value no, 
and where zo = 1 + B^^'> + B^^\ In Fig. O we have 
shown the density and fluctuations as a function of the 
chemical potential. With increasing temperatures, the 
step-like density profile characteristic of a Mott state be- 
comes increasingly smoother, reflecting the fact that the 
free energy cost to create supplementary particles or holes 
vanishes near the edges of a Mott plateau. A Mott-like 
region survives up to T* ~ 0.2 U/kB (see also 0]), which 
can be seen as a melting temperature for the Mott phase. 
We have also compared the exact expressions to the PHA, 
without noticeable difference up to T* . 

With the properties of the zero-tunneling model clari- 
fied, we reintroduce tunneling in a second step. At zero 
temperature, when exploring a Mott domain at constant 
J and varying chemical potential, the Mott state remains 
stable as long as its elementary (quasi-particle and hole) 
excitations are. At the upper or lower Mott lobe bound- 
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FIG. 2: Density profile and fluctuations for an array of iso- 
lated wells (J = 0), shown for various temperatures. 



aries, the chemical potential /i*^^^ becomes equal to the 
excitations energies £'^'^0 '^'^'^ , thus favoring the prolifera- 
tion of particle or hole excitations to reduce the free en- 
ergy [l9| . To investigate how this behavior changes with 
increasing temperature, we have calculated the quasipar- 
ticles and holes dispersion relations at finite temperature 
using the PHA supplemented by a random-phase approx- 
imation as in 0, Sim, [111. We find 



(3) 



with Ak = ^C'^Jl + 2UDJ]^ + U'^, Jk = -2zJxk and 
Xk = 7 Si=i '^^'^ finite temperature disper- 
sion relations are formally similar to the zero temper- 
ature expressions [2^ except for the factors C/D = 
(no + l)^(^P) T noA^'^^\ with A^^p/^'^) = (1 - S(±))/zo. 

In the following, we mostly work for simplicity in the 
limit U 3> zJ, where the dispersion relations reduce 



to E, 



(qp/qh) 



(qp/qh) 




=p j(qp/qh)^^^ ^ith an eflFective 
tunneling energy J^'^p' = zJ{no + 1)^('5p) for particles 
and jf'^'') = zJno^''^''-' for holes. We concentrate first 
on the upper boundary and try to find the chemical 



potential for which E, 



(qp) 

k=0 



(+) 



Introducing a vari- 



able s ~ {Uno — /.j^+')/A:bT, we find that, up to terms 
~ e^^^ , s solves the equation s — atanh (|) = 0, with 
a ~ zJirio -\- l)/kBT. Graphic inspection shows that 
there is no solution for a < 2, corresponding to a critical 
temperature ksTc = zJ(no + l)/2 above which there is 
no SF region left due to thermal depletion. This argu- 
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ment holds far from the zero-temperature Mott transi- 
tion, where Eqs. ([3]) are vaHd. To obtain the complete 
phase boundaries, we have solved numerically the equa- 



tions 



(qp/qh) 
k=0 



/i^*) for any zJ /U . The result is shown 
in Fig. [U where three distinct regions emerge. For 
zJ <^ 2kBT / (jiQ + 1), one goes continuously from one 
insulating "lobe" to the other through a region in the 
normal phase. For 2kBT/{nQ -I- 1) -C zJ <C zJc, well- 
defined Mott lobes survive, as in the ground state, and 
above the critical Jc a phase transition takes place to a SF 
phase, albeit with reduced SF fraction. Above T ^ T* , 
we expect the system to be essentially normal for J < U. 
Note that when Mott regions with different filling factors 
due to an external trap potential, the regimes described 
above can coexist. 

To characterize the thermodynamics for x < T* 
more quantitatively, we calculate the free energy as 
F = Fo + AF('iP) + AF^i^\ The first term Fq = 
— fcsTln(l -I- + gives the free energy corre- 

sponding to the thermal activation of local defects in 
each well for zero tunneling. The remaining contribu- 



tions AF(qp) 

^^(qh) ^ 



Fi 



BE 



F, 



(qp) 



Fbe [U no - fi] and 



- Fbe [/i - U{no - 1)] de- 



scribe instead dilute gases of quasi-particles and quasi- 
holes mobile through the lattice. Here, the functional 
-FbeN = X^k (1 ^ cxp(/3a;)) gives the free energy of an 
ideal Bose gas. For C/ 3> J, we obtain for instance 



E 



V 



V 



1(4,) 



and a similar expression for AF^'ih). Here Iq is a Bessel 
function, and 2? the dimensionality. 

Thermodynamic quantities follow from derivatives of 
the free energy. Contributions from local activation, 
denoted by the index "o" have been calculated above. 
Terms coming from particle and hole contributions 
are readily obtained. For instance, the particlc/holc 
corrections to the average density read An^'^P/'^''^ = 
— 9AF*^'^P/''^-'/9/i, whereas the density fluctuations can 
be characterized using the Hellmann-Feynman theorem 
[2^, which to leading order in e^^^ , reduces to Var(n) w 
Var(n)o + An('iP)(l - An^^P)) - ArT('5h)(i + An^'i^'>). 

In the low-temperature limit, T <^ Tc, the contribu- 
tion of local activation terms is negligible, and tunnel- 
ing corrections to the thermodynamic functions domi- 
nate. However, using the asymptotic form for Iq we 
find for instance AF^ip^ /NsksT = [kBT/iTrJ{no + 

l)]'/'ff5/2(e''^'^"^*^')), where = de- 

notes a Bose function. The prefactor implies that cor- 
rections to the zero-temperature behavior are highly 
suppressed. Hence this regime correspond to a quan- 
tum Mott insulator, where the properties of the zero- 
temperature system are essentially preserved up to small 



thermal corrections. Increasing the temperature, correc- 
tions from local activation and from tunneling terms be- 
come progressively comparable until the high tempera- 
ture regime, T ^ Tc, is reached. In this regime, using 
the Taylor expansion Jo (a;) ~ 1 -f for the Bessel 

function, we find 



NskeT 



1 

'4^ 



J{no + l )y F(+) 
keT 



(5) 



Note that this expression also applies to the normal phase 
in-between two Mott-like regions. Here the right hand 
side has a small term oc (Jno/ZcsF)^ <C 1, times a factor 
comparable to the J = result. This means that in this 
temperature regime, the model with J = alone already 
gives a good description of the thermodynamics. 

As explained above, the results can be directly applied 
to the harmonically trapped case in the LDA. We ap- 
ply now the calculation to estimate the temperatures 
achieved in current experiments. Typically, a Bose- 
Einstein condensate is first produced in a magnetic trap 
with frequency uji, at an initial temperature T^. We cal- 
culate its entropy using the Popov approach as exposed 
in [2^ 1 using the local density approximation and for typ- 
ical experimental values = 27r x 20 Hz and N = 2x 10^ 
atoms. Then, the cloud is slowly transferred in the opti- 
cal lattice of depth Vq. This also changes the frequency 
LUf of the external harmonic potential according to 



8^0 



(6) 



In Eq. © the second term comes from the trapping force 
due to the Gaussian shape of the laser beams forming 
the lattice, with w the lattice laser waist (1/e^ radius), 
typically w 150 fim. Assuming that the transfer into 
the lattice is adiabatic, both entropy and atom number 
are conserved. In the final state, assumed to be deep 
in the MI regime, the number density n = —df/d^ and 
entropy density s = —df/dT are calculated from the 
free energy density / [2^. Those densities are then inte- 
grated over space to fix the total atom number and en- 
tropy, which determines the final chemical potential and 
final temperature Tf. We plot in Fig. [3^ the result of a 
self-consistent calculation solving for Tf for a fixed atom 
number N = 2 x 10^. One sees that for small initial tem- 
perature Ti < 0.4 Tco (condensed fraction larger than 75 
%) one ends well within the Mott region. Only by start- 
ing from a rather large thermal fractions does the final 
temperature exceed 0.2 U. Reaching the quantum insu- 
lator regime is more challenging, but still within reach of 
current experimental possibilities {Ti < 0.2 Tco). 

For this calculation, the role of the trapping poten- 
tial is essential. Similar calculations for homogeneous 
systems with integer filling would yield a much higher 
temperature [l^, since near the center of the Mott lobe 
the entropy is exponentially suppressed by the finite in- 
teraction energy U. In the trapped case, the entropy in 
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fact concentrates around the lobe boundaries, where the 
generation of excitations is easiest. As long as T < T*, 
the peak value of the entropy density is s^ax ~ log(2)- 
The thickness 6R of the layer where thermal fluctuations 
are important can be estimated from moj^RoSR ~ k^T, 
where Rq = y^2U / muj'^cP is a typical spatial scale 
for the shell structure. Since the total entropy and 
atom number scale respectively as S'/fce ^ 47ri?Q5i? and 
N ~ {AT:/3){Ro/d)^^^, we find that in the trappe case, 
the entropy per particle rises linearly with temperature, 
S/NkB ~ SfceT/f/, instead of exponentially (- e-^l^^"^) 
as in the uniform case (see also [l2| for a similar analysis). 

An important consequence of Eq. ([5]) is that the trap 
frequency increases significantly when the lattice depth 
is raised. This leads to spatial compression (and the for- 
mation of Mott plateaux with higher filling) , but also to 
adiabatic heating of the cloud, as in a conventional har- 
monic trap. In the second calculation, this effect was 
investigated [2^. In Fig. [3Jd, the result is plotted for a 
fixed initial temperature T,; = 0.3 Tco, but varying final 
trap frequencies. It is clear that lower and lower temper- 
atures are achieved when the external trap frequency is 
reduced. In particular, one can achieve T < 0.01 U if 
the transfer is done at constant trap frequency. Exper- 
imentally, this could be realized using bluc-dctuncd lat- 
tice beams, or additional bluc-dctuncd (repulsive) lasers 
to compensate for the trapping force of the laser beams. 
Note finally the bump in the curve, that indicates that 
the transition from a unity-filled MI to a two-plateaux 
distribution with cither single or double occupancy. 

In conclusion, we have discussed the thermodynamics 
of ultracold bosons in optical lattices. We have identified 
two characteristic temperatures, T ~ Tc, above which the 
SF regions disappear, and a melting temperature T* « 
0.2 U. The thermodynamic function for the Mott phases 
arc calculated explicitly. We have used the results to 
estimate the temperature reached in current experiments, 
and found that they easily reach the thermal insulator 
regime T < T* , and possibly the quantum region T < 
Tc- We suggest an adiabatic decompression scheme which 
potentially allows to reach much lower temperatures in a 
system with unity filling. 
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FIG. 3: Final temperature achieved by adiabatically trans- 
ferring a condensate into the optical lattice potential (depth 
Vb = 20 En)- The transfer is done in (a): for a fixed final trap 
frequency lj/ = 27r x 70 Hz varying the initial temperature 
Ti, and in (b): for a fixed initial temperature Ti = 0.3 Tco 
but varying the final trap frequency ujf. The solid line in- 
dicates the melting temperature T*. Above the dashed and 
dotted lines, the superfluid layers corresponding to no = 2 and 
no = 1, respectively, turn normal. For both plots, the initial 
frequency coi = 2n x 20 Hz and the atom number N = 2x 10^ 
are fixed. 
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